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Abstract 

The aim of this note is to prove a resuh of efTective stabihty for a non-autonomous perturbation of an integrable 
Hamiltonian system, provided that the perturbation depends slowly on time. Then we use this result to clarify 
and extend a stability result of Giorgilli and Zehnder for a mechanical system with an arbitrary time-dependent 
potential. To cite this article: 

Resume 

Stabilite effective pour des Hamiltoniens presque integrables lentement non-autonomes et applica- 
tion. Le but de cette note est de demontrer un resultat de stabilite effective pour une perturbation non-autonome 
d'un systeme hamiltonien integrable, sous la condition que la perturbation depende lentement du temps. Nous 
utilisons ensuite ce resultat pour clarifier et generaliser un resultat de stabilite de Giorgilli et Zehnder pour des 
systemes mecaniques dont le potentiel depend arbitrairement du temps. Pour citer cet article : 



1. Introduction 

Let n e N, n > 2, T" = M"/Z" and consider the Hamiltonian system defined by i? : T" x K" x M ^ M, 

H{6,I,t)^h{I)+ef{e,I,t), (0,/,^) = (^^l,...,(?„,/l,...,/„,^)eT"xM"x]R, e>0. (1) 

Nekhoroshev proved ([1]) that whenever h is steep (see §2 for a definition), f{6,I,t) — f{0.1) is time- 
independent and H is real-analytic, there exist positive constants eo, ci, C2, C3, a, 6 such that for e < £0, 
for all solutions {9{t),I{t)) we have the following stability estimates 

- /(O) I = max |/,(t)- 7,(0)1 <Cle^ |i| < C2 exp(c3e-'^). (2) 

l<i<n 
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In the particular case where h is (strictly uniformly) convex or quasi-convex, following a work of Locliak 
([4]) it was proved ([2], [3]), using preservation of energy arguments, that one can choose a = b — (2n)^^ 
in (2), and that these valiics arc close to optimal. In the general steep case, however, there are still no 
realistic values for these stability exponents a and b. The purpose of this note is to discuss to which extent 
stability estimates similar to (2) hold true if the perturbation is allowed to depend on time. 

Assume first that / depends periodically on time, that is f{9, 1, t) = f{9, 1,t + l) in (1). Removing the 
time-dependence by adding an extra degree of freedom, the Hamiltonian is equivalent to 

H{0, if, I, J) = h{I, J) + ef{0, <p, I), {0, ^ = i, /, J) G T" X T X R" X R, h{I, J) = h{I) + J. 

It turns out that if h is convex, then h is quasi-convex and so (2) holds true with a = b = (2(n-|- In 
general, it is possible for h to be steep in which case (2) is satisfied, but it is not clear how to formulate 
a condition on h (and not on h) to ensure that (2) holds true. 

Now assume that / depends quasi-pcriodically on time, that is f{0,I,t) = f{0,I,tLu) in (1) for some 
vector oj G which we can assume to be non resonant (A; • w 7^ for any non-zero k G Z™) and 
/ : T" X R" X T™ — >■ R. As before, the time-dependence can be removed by adding m degrees of freedom 
and we are led to consider H{0, ip, I, J) = h{I, J) + sf{0, (p, I) but this time 

(61, <p = tui, I, J) G T" X X R" X R", h{I, J) = h{I) + w • J. 

It was conjectured by Chirikov ([5]), and then again by Lochak ([6]), that if h is convex and ui satisfies 
a Diophantine condition of exponent r > m — 1 (there exists a constant 7 > such that \k.uj\ > j\k\~'^ 
for any non-zero k G Z"^), then the estimates (2) hold true and moreover we can choose a = b = 
(2(n + 1 + t))~^ . If m = 1, then r = and we are in the periodic case so the conjecture is true. However, 
if m > 1, ft, cannot be steep and the problem is still completely open. Even though the conjecture is 
sometimes considered as granted (for instance in [7]), there is still no proof. Needless to say that the 
situation in the general case (without the convexity assumption on h) is even more complicated. 

In a different direction, Giorgilli and Zehnder ([8]) considered the following time-dependent Hamiltonian 

G{0,I,t) = h2{I) + V{0,t), (^,7,i) gT" xR" xR, h2{I) = + ■ ■ ■ + 1^, 

and proved the following Nekhoroshev type result: if G is real-analytic and V is uniformly bounded, then 
for R sufficiently large, if Iq belongs to the ball Bn of radius R centered at the origin, then I{t) G B2R 
for \t\ < C2 exp(c3_R'^) for some positive constants 02,03 and d. Even though such a system is clearly 
not of the form (1), the fact that no restriction on the time-dependence is imposed in their result has 
lead to several confusions. In [8], the authors themselves assert that "extra work is needed because the 
time-dependence is not assumed to be periodic or quasi-periodic" . Even more surprising, one can read (in 
[9] for instance) that this result implies that the estimates (2) hold true for (1) without any restriction 
on the time-dependence. Concerning the latter assertion, it is simply wrong and it seems very unlikely to 
have non trivial stability estimates for (1) with an arbitrary time-dependence. As for the former assertion, 
it is not difficult to see that the system considered in [8] can be given the form (1), with a perturbation 
depending slowly on time, in their example f{9,I,t) = f{0, 1,e^^^t). We will show in §2 that for such 
Hamiltonian systems depending slowly on time, essentially classical techniques can be used to prove 
that (2) hold true, and that the non-periodicity or non-quasi-periodicity of time plays absolutely no role 
(as a matter of fact, we already explained that for a periodic or quasi-periodic time dependence which is 
not slow, basic questions are still open). Then, in §3, we will use this result to derive, in a simpler way, a 
more general statement than the one contained in [8]. 
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2. A stability result 



For a given p > 0, a function h G C'^{Bp) is said to be steep if for any k € {1, • ■ • ,n}, there exist 
positive constants Pk,Ck,Sk such that for any afEne subspace Afe of dimension k intersecting Bp, and 
any continuous curve 7 : [0, 1] ^ Afe n B with |7(0) — 7(1)1 = S < Sk, there exists e [0, 1] such that 
ItW ~ 7(0)1 < ^ * ^ [0,t*] and |nAj,(V/i(7(i*)))| > 0^5'^'', where IIa;, is the projection onto A^, 

the direction of Afe. Then, given r, s > 0, let us define the complex domain 

Vr,s = {{e,I,t) e (C"/Z") X C" X C I |(Im(^i), . . . ,Im(^„))| < s, \lm{t)\ < s, d{I,Bp) < r), 

and for a fixed constant 1/2 < c < 1, we consider H{9, /, t) = h{I) + ef{9, /, eH) defined on P^.s and real- 
analytic (that is H is analytic and real- valued for real arguments). Finally, we assume that there exists a 
positive constant M such that the operator norm |V^/i(/)| < M for any / e B^, and that \f{9,I,t)\ < 1 
for any (9,I,t) G I?r,s- 

Theorem 2.1 Under the previous assumptions, there exist positive constants eo,ci, 02,03, that depend 
on n, p,pk,Ck, Sk,r, s, M , and positive constants a,b that depend only on n,pk, such that if e < Eq, for 
all solutions {9{t),I{t)) of the Hamiltonian system defined by H, if I{Q) € Bp/2, then \I{t) — /(0)| < cie'' 
for all \t\ < C2 exp(c3e~"). 

Let us explain the proof, which follows from the arguments exposed in [10] or [11], up to some technical 
points we shall detail now. First we remove the time-dependence: we let x = eH and we introduce a 
variable y canonically conjugated to x, so that the Hamiltonian is equivalent to 

H{9, 1, X, y) = h{I) + e'y + ef{0, 1, x) = h{I) + f{0, 1, x, y), {9, 1, x, y) e (3) 

where V^^s = T^r.s x {?/ G C | |Im(7;)| < s}. The fact that the dependence on time is slow allows us 
to keep the integrable part fixed when removing the time-dependence, as one can consider that the 
extra degree of freedom only affects the perturbation. The new perturbation / depends on parameters 
or "degenerate" variables x and y (degenerate since they are not present in the integrable part), and 
such systems were already considered by Nekhoroshev ([1]). However, a difficulty arise: for subsequent 
arguments, it is important for the (real part of the) variable y to be unbounded, which is indeed the case 
by our definition of but on this extended domain / is unbounded and this implies that H in (3) 
is not a perturbation of h. Yet the Hamiltonian vector field Xjj can be considered as a perturbation 
of Xh, as Xj = {dif, -def, dyf, -d^f) = {edjf, -edef, e", -dtf), and so X^ is bounded (by a Cauchy 
estimate) on the domain 'Dr/2,s/2 for instance. As a consequence, even when h is convex one cannot use 
preservation of energy arguments as it is the case in [2] , [3] , [4] , and in general one has use a pcrtiubation 
theory that deals only with vector fields: the proofs in [10] and [11] accommodate both requirements. 
Now the analytic part of [10] and [11] goes exactly the same way for (3) by simply considering .7: and y as 
"dummy" variables: given an integer parameter m > 1 which will be determined by the geometric part in 
terms of e, on suitable domains resonant normal forms with a remainder of size bounded by a constant 
times e'^e"'" are constructed (note that the size of the perturbation X j is of order e"^ and c < 1, but 
its "effective" size is of order e and so m will be determined in terms of e and not e'^; e'^ just enters the 
pre- factor in the exponential and will not alter the radius of confinement as we always have 6 < 1/2 
whereas c > 1/2.). The geometric part of [10] and [11] also goes exactly the same way since the time of 
escape (of the domain) of the degenerate variables x and y is infinite (as the domain is unbounded in these 
directions), m is eventually chosen proportional to e~", the radius of confinement is chosen proportional 
to £^ and the stability time is bounded by a constant times e"™. 

Now let us add two remarks on the statement of Theorem 2.1. First, the exponents a and b are the 
same as in (2) when the perturbation is time-independent. It is reasonable to expect that if h is convex, 
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then a = b = (2n)~^ in Theorem 2.1, but we already explained that we cannot use preservation of energy 
arguments and so we cannot reach these values: the problem actually reduces to the problem of finding 
realistic values of a and b in the general steep case, which is still open. Then, using [10] and [11], the 
statement of Theorem 2.1 can be generalized in two ways: using [10] the statement holds true for the 
much wider class of Diophantine steep functions introduced by Niederman (which is a prevalent class of 
functions), using [11] the statement holds true for a-Gevrey Hamiltonians for a > 1 (with exp(c3e~'') 
replaced by exp(c3e~" recall that 1-Gevrey is real-analytic) and for C'^ Hamiltonians, k > n + 1 
(with exp(c3e~") replaced by cse"'^ with k* the largest integer I > 1 such that k> ln+ 1). 



3. An application 

Now we come back to the problem studied in [8], and more generally we consider, for an integer p > 2, 

Gie,I,t) = hpiI) + Vie,t), (e,7,i) gT" xR" xM, /ip(7) = 7f + • • • + 7^. 

The case p — 2 corresponds to [8] and h2 is convex, for p > 3 the function hp is not convex but it is 
steep with p/^ = p — 1 and Ck = (5^ = 1 for all 1 < fc < n. The function V is assumed to be real-analytic, 
defined on Vs = {{e,t) e (C"/Z") x C | |(Im(ei), . . . ,Im(e„)| < s, \Im{t)\ < s}, and it is assumed that 

\V{9,t)\ < 1 for all (9,1) G Vs. 

Theorem 3.1 Under the previous assumptions, there exist positive constants 7?o,ci,C2,C3 that depend 
on n,p,s, and positive constants a',b' that depend only on n,p, such that if R > Ro, for all solutions 
{0{t),I{t)) of the Hamiltonian system defined by G, if I{0) G Br, then \I{t) — 7(0)| < ciR^~^ for all 

\t\ < C2i?i"Pcxp(c3i?"'). 

The proof is a direct application of Theorem 2.1. Indeed, for i? > consider the scalings 

7 = 7?7', 9 = e', G = RPG', t = R'^'^t' . 

Then the Hamiltonian G{6, 1, t), for {$, 7, f) e Dg x B2R, is equivalent to the Hamiltonian G'{6', I', t'), for 
ie',r,t') eVsX B2, where G'{e',r,t') = hpil') + R-PV{e',R^-Pt'). Hence we can apply Theorem 2.1 
to the Hamiltonian G', with e = R~p, c = {p — l)p~^, p = 2 and M which depends only on p: there 
exist positive constants eo,ci,C2,C3, that depend on n,p,s, and positive constants a,b that depend only 
on n,p, such that if £ < Eq, for all solutions {8' (t'), I' {t')) of the Hamiltonian system defined by G", if 
7'(0) e Bi, then {!' {t') - 1' {0)\ < as'' for ah \t'\ < C2 exp(c3e"''). Recalling that e = R-p, this means that 
if 7? > 7?o = e/"\ for all 7'(0) e Bi we have |7'(f )-7'(0)| < ciR-'>' for all \t'\ < C2 cxp(c3i?"') for b' = pb 
and a' = pa. Now scaling back to the original variables, for all 7(0) G Br, we have \I{t) — 7(0) | < CiR^~^ 
for all \t\ < C2i?i-Pexp(c3i?"'). 

Now let us add some comments on the statement of Theorem 3.1. The estimate \I{t) — 7(0) | < ciR^~^ 
is stronger than |7(t) — 7(0) | < R (as in the argument above, the estimate |7'(t') — 7'(0)| < Cie^ is 
stronger than \I'{t') — 7'(0)| < 1) and hence it is stronger than I{t) € B2R if 7(0) € Br. Moreover, we 
have \t\ < C2exp(c3i?" ) < C2-R^^p exp(c3i?° ) by restricting C3 to a smaller value C3 and enlarging Rq if 
necessary. So even for the convex case p ~ 2 our statement is more accurate than the statement in [8]. 
In fact, for p = 2, we already explained that wo believe we can choose a = b = (2n)^^, in which case the 

statement of Theorem 3.1 would read \I{t) — 7(0)| < ciR^~" ^ for all \t\ < C2R^^ cxp (^c^R" , which 
would be in perfect agreement with the much simpler autonomous case V{9,t) = V{6) described in [3]. 
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